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electric field in its vicinity. This leads to vacuum instability and to the for- 
mation of a charged vacuum condensate which screens the electric charge of 
the string. We analyze the structure of this condensate using the Thomas- 
Fermi method. 
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I. INTRODUCTION 

Cosmic strings are linear defects that could be formed at a phase transition in the early 
universe. (For a review see []]].) Witten has shown that strings predicted in some grand 
unified models behave as superconducting wires. Such strings moving through magnetized 
cosmic plasmas can develop large currents and can give rise to a variety of astrophysical 
effects. In particular, they have been suggested as possible sources of ultrahigh energy 
cosmic rays ||. 

Currents developed by oscillating strings in a magnetic field are not homogeneous along 
the strings because different portions of the string cross the magnetic field lines in different 
directions. This results in charge accumulation, and portions of the string can develop a 
charge per unit length A comparable to the current, A ~ J . (Here and below we use units 
in which h — c — 1.) The electric field near the string is given by 

B=%. (1) 

r 

It can become extremely strong in the immediate vicinity of the string, and then quantum 
effects, such as vaccum polarization and pair production, must be taken into account. One 
can expect that the created particles of charge opposite to that of the string will accumulate 
in bound states and form a condensate screening the electric field near the string to below the 
critical value. It was noted earlier that such screening would lead to a drastic modification 
of string electrodynamics f|] and would have a significant effect on the propagation of high- 
energy particles emitted from the charged portions of the string ||. The purpose of the 
present paper is to give a quantitative description of the screening condensate near a charged 
superconducting string. 

The superconducting current in the strings is carried by charged particles which acquire a 
mass M at the string-forming phase transition but remain massless inside the strings. These 
massless charge carriers move along the strings at the speed of light. The string current is 
bounded by the critical value, 



J c ~ eM, (2) 

at which the characteristic energy of the charge carriers becomes comparable to M, so that 
they have enough energy to jump out of the string. The mass M is model-dependent but is 
limited by the string symmetry breaking scale i], M < 77. 

In a cosmological setting, the string charges and currents vary on astronomical time and 
length scales, and for our purposes we can regard them as constant. For a string segment 
with J < A, we can always find a Lorentz frame where J = 0. The electric field close to the 
string will then be well approximated by that of an infinite straight string. We shall consider, 
therefore, an infinite straight string with a constant charge per unit length A and vanishing 
current, J = 0. For sufficiently large A, charged particles (for definiteness electrons) have 
bound states localized near the string with negative energies smaller than — m, where m 
is the electron mass. The vacuum then becomes unstable with respect to production of 
electron-positron pairs. For a positively charged string, positrons are repelled away, while 
electrons form a vacuum condensate surrounding the string. We shall determine the electric 
field and the charge distribution in this condensate using the Thomas-Fermi method |5| 
in which the condensate is approximately treated as an ideal gas obeying the Fermi-Dirac 
statistics. 

In the next section we shall review the derivation of the relativistic Thomas- Fermi equa- 
tion and specify the boundary conditions appropriate for the case of cylindrical symmetry. 
Approximate analytic solutions of this equation are given in Sec. Ill, and its numerical solu- 
tions are presented in Sec. IV. The conclusions of the paper are summarized and discussed 
in Sec. V. 

II. THOMAS-FERMI EQUATION 

The density of electrons in a degenerate Fermi gas is related to the Fermi momentum pp 
by 



n e (r) 



Pf 
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(3) 



The relativistic relation between the Fermi energy ep and Fermi momentum is 
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where m is the electron mass, — e is its charge, V(r) = —e(p(r) and (p(r) is the self-consistent 
electrostatic potential for an electron, taking into account both the field of the string and 
the average field produced by other electrons of the condensate. The condensate is formed 
of electrons occupying quantum states in the negative energy continuum, e < —m. We 
therefore set the fermi energy to be €f = —m. The electron density (|3]) is then given by 



n e (r) 



1 
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(5) 



Introducing the total charge density pr which is composed of the electron charge and 
external string charge, 



Pt = p a — en e 



(6) 



and using the Poisson equation 



AV(r) = 47rep T (r), 



(7) 



we find a self-consistent non-linear differential equation 



AV(r) = -47re 
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(y 2 {r) + 2mV{r)) - p s (r) 



(8) 



This equation has been used in [|M to study the electron condensate around supercharged 
nuclei. In the case of a string, the problem has cylindrical symmetry and 



AV(r) = V"(r) + -V'(r). 
r 



(9) 



We shall approximate the string charge distribuition as a uniform distribuition in a 
cylinder of radius S, 



p s (r) = p 9(5-r). (10) 

The linear charge density of the string is given by A = 7r<5 2 po- The charge carriers are 
typically concentrated in a tube of radius r ~ M _1 ; hence, we should have 5 ~ M~ 1 . 

It is easily seen from Eq. (|J) that the density of electrons, n e (r), is different from zero 
only in the region of space where V(r) < —2m. Therefore, the condensate has a finite radius 
r = R c . For r > R c , the solution of (§) is just the usual logarithmic potential of a linear 
charge, 

V(r) = 2e\ \n^- (11) 

-ft* 

Here, Ao is the total charge per unit length of string, including both the charge carriers in 
the core and the condensate, and R* is the cutoff radius indicating the distance at which the 
approximation of an infinite straight string breaks down. R* is given by the smallest of the 
following three length scales: (i) the typical distance between the strings in a cosmic string 
network, (ii) the characteristic curvature radius of string, (iii) the typical wavelength of the 
current-charge oscillations along the string. 

The boundary condition for Eq. (^) at r = is 

HO) = 0, (12) 

while at r = R c we have 

V(R C ) = -2m, V W = n^km - (13) 

Note that we have three rather than two boundary conditions, as a second-order differential 
equation would normally require. The third condition is needed to determine the condensate 
radius R c 

We expect R c to be microscopic, while i?* will typically be astrophysically large. Hence, 
the logarithm in Eq.(|13D is ln(i?*/i? c ) ~ 10 2 . In numerical calculations below we choose i?* 
so that ln(i?*/i? c ) ~ 30; our results are not sensitive to this choice. 



The Thomas-Fermi approximation is adequate when the characteristic scale of variation 
of the condensate density n e (r) is large compared to the electron wavelength l/p(r). The 
corresponding condition is 

1 



d 

dr 



< 1. (14) 



p(r) 

We shall see that this condition is satisfied in most of the condensate region < r < R c , 
provided that the charge density A is sufficiently large. 

III. ANALYTIC APPROXIMATIONS 

The Thomas-Fermi equation (|]) can be solved analytically in the limit when the mag- 
nitude of the potential V(r) is large, |V(r)| ^> 2m. We can then neglect 2mV(r) compared 
to V 2 (r), and outside the string core Eq.(H) reduces to 

V n (r) + -V , (r) = ~\V(r)\ 8 . (15) 

T 07T 

This has a solution 

V{r) = -C/r (16) 

with 

C = (37r/4e 2 ) 1/2 « 18. (17) 

The corresponding electric field is 

E(r) = C/er 2 . (18) 



We note that the solutions ( JIB]) and (0) do not depend on the string charge density A. 
As r decreases, the electric field ( Pf ) grows faster than that of the vacuum solution ([I]). It 
cannot, therefore, be extended all the way to the string but has to be matched with Eq. ([!]) 
at some radius R s below which the vacuum solution takes over. The matching radius at 
which the two electric fields become comparable is 
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R s ~ C/eX ~ 200A" 1 . (19) 

We shall call it the screening radius. For r <C R s , the screening is unimportant and the 
electric field is given by Eq. (jT]) . The screening radius is always large compared to the string 
thickness 5 ~ M _1 , provided that A is smaller than the critical value (Q), 

A < eM. (20) 

The potential corresponding to the vacuum solution ([!]) at 5 < r <C R s is 

V(r) = -2eX[ln(R s /r) + B], (21) 

where B ~ 1 is a numerical constant. The potential at the string core is thus 

V(0) « -2eAln(i? s /5). (22) 

The condition |V"(r)| ^> m implies r <C C/m, and thus the solution (|T8|) is valid in 
the range C/eX Cr <C C/m. This range exists only if A is sufficiently large, A ^> m/e. 
Combined with the condition ( p0|) this implies M 3> m/e 2 . In models of astrophysical 
interest, the charge carrier mass M is very large (so that the strings can develop large 
currents and charges), and this condition is satisfied with a large margin. 

At r ~ C/m, the potential V(r) becomes comparable to — m signalling that we are close 
to the condensate boundary [see Eq. (JT~5|)] . Hence, we can estimate the condensate radius as 

R c ~ C/m. (23) 



The condition of validity of the Thomas- Fermi approximation (14]), when applied to the 



solution (|16D, gives C ^> 1. This is satisfied with a reasonable accuracy [see Eq. ([P7D1. 



IV. NUMERICAL CALCULATION 



We obtained numerical solutions to the Thomas- Fermi equation for Eq. (|S[) for V(r) 
with the boundary conditions fll2|) and (|13|) using the relaxation method. The resulting 



electric field is plotted in Fig. 1, together with the analytic approximations ([!]) and (|T8|). 
The agreement between the analytic and numerical solutions is excellent in the appropriate 
ranges of the radius r. 

We have verified that the shape of V(r) outside the string core is not sensitive to the 
value of the core radius 5. In particular, the condensate radius R c approaches a constant 
value independent of S (see Fig. |2|). This is very fortunate, since a realistic value of the 
core radius would be too small to resolve in our calculations. Figure || suggests that it is 
sufficient to choose 5 <C m _1 . We used 5 = 10 _3 m _1 in most of the calculations described 
below. 

The condensate radius R c is plotted in Fig. ^ as a function of the linear charge density 
of the string, A. We see that at large A, R c approches a constant value, 

R c = 82m" 1 , (24) 

in agreement with Eq. (|23|). The screened linear charge density of the string Ao, which 
determines the electric field outside R c , is also found to be independent of A: 

A « 5.34em. (25) 

The electric field at the condesate boundary is 

E = 2\ /R c « 10 _1 em 2 . (26) 

Note that this is considerably smaller than the critical field, E c = m 2 /e, which signals the 
onset of intensive pair production [|IIJ]. In our case, E Q ph 10~ 3 E c . We shall return to this 
point later in Sec. V. 

The effective linear charge density A e //(r) inside the condensate can be found as 

Wrl.J,^*''— (27) 

where we have used Eq. (|7|). As r grows, A e //(r) decreases and we can define the effective 
screening radius R s as the radius at which half of the string charge is screened, 
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\ eff (R s ) = A/2. (28) 

At the boundary of the condensate we must have X e ff(R c ) = Xq. The screening radius R s 
is plotted in Fig. 4 for several values of A. We see that, although the condensate radius R c 
is independent of A, the screening radius gets smaller rather rapidly as A is increased. A 
numerical fit to the data in the Fig. 4 gives 

R s = 80A"\ (29) 

in agreement with the order-of- magnitude estimate (|19|). 



V. SUMMARY AND DISCUSSION 

We have found that a superconducting cosmic string having a sufficiently large charge 
per unit length, A ^> m/e, is surrounded by an electron condensate of radius R c ~ lOOm -1 , 
where m is the electron mass. In the immediate vicinity of the string, the effect of the 
condensate is unimportant and the electric field is given by the vacuum solution, E m 2X/r. 
Screening due to the condensate becomes significant at r ~ R s ~ 100A -1 , and for R s <^ r <^ 
R c the electric field has the form E « (37r) 1//2 /2e 2 r 2 . Outside the condensate, at r > R c , the 
field is given by E ~ em/r « 10 -2 em 2 (R c /r). 

As we already mentioned, the electric field at the condensate boundary is well below the 
critical field, E ~ 10~ 3 E C , where E c = m 2 /e. The rate of pair production per unit volume 
in a homogeneous electric field is [[IT]] 

dN/dVdt « (eE/ir) 2 exp(-7rE c /E), (30) 

which indicates that the outer parts of the condensate where E <^C E c will be filled up very 
slowly. The characteristic time of pair production, r ~ (eE) -1 ' 2 exp(%E c / E), is greater 
than the age of the universe for E ^$ 4 x 10 2 -E C . For astrophysical strings, we expect the 
condensate radius to be given by the distance from the string at which such values of the 
electric field are reached. From Eq. (|16|) we find 



R c ~ —m~ l « 0.2m- 1 . (31) 

20 v ; 

As the charge density of the string A is increased, the potential near the string core 
becomes more and more negative. As a result, particles more massive than electrons de- 
velop condensates. From Eq. (p2|), particles of mass \i develop states with e < — fi at 
A ~ fi/eln(R s /5) ~ \x. The condensates of different particle species will have the form of 
coaxial cylinders, with condensates of more massive particles being closer to the string. 

Finally, we would like to mention some open questions. In this paper we studied vacuum 
condensation of fermions. Charged Bose particles, such as Higgs and gauge bosons will also 
form vacuum condensates, and the properties of these bosonic condensates may differ from 
the fermionic case. Another important problem is the nature of modifications introduced 
by vacuum screening in string electrodynamics and in the propagation of charged particles 
emitted by the strings. We hope to return to some of these issues in future publications. 

ACKNOWLEDGMENTS 

J. R.S.N, is grateful to the Institute of Cosmology, Tufts University, for hospitality. The 
work of J. R.S.N, was supported in part by funds provided by Conselho Nacional Desenvolvi- 
mento Cientifico e Tecnologico, CNPq, Brazil. The work of I.C. and A.V. was supported in 
part by the National Science Foundation. 



10 



REFERENCES 

[1] A. Vilenkin and E.P.S. Shellard Cosmic Strings and other Topological Defects ( Cam- 
bridge University Press, Cambridge, England, 1994). 



[2] 
[3] 

[4] 
[5] 
[6] 
[7] 



[9 
[10 



E. Witten, Nucl. Phys. B 249, 557 (1985). 

C.T. Hill, D.N. Schramm and T.P Walker, Phys. Rev. D 36, 1007 (1987). 

M. Aryal, T. Vachaspati and A. Vilenkin, Phys. Lett. 194 B, 25 (1987). 

V. Berezinsky and H.R. Rubinstein, Nucl. Phys. B 323, 95 (1989). 

L. H. Thomas, Proc. Camb. Phil. 23, 542, 1927; E. Fermi, Z. Phys. 48, 73 (1928). 

M. S. Vallarta and N. Rosen, Phys. Rev. 41, 708 (1932). 

A. B. Migdal, V. S. Popov and D. N. Voskresenskii, Sov. Phys. JETP 45, 436 (1977). 

B. Miiller and J. Rafelski, Phys. Rev. Lett. 34, 349 (1975). 
J. Schwinger, Phys. Rev. 82, 664 (1951). 



11 



FIGURES 



20 



15 



LU 
LU 



-5 



-10 



-15 







I I I 


I 




-y/ 


8 






- 


- 










- 


- 




R \ 

S \ 




- 








\\ 




- 






\\ 




- 








\ \ 












v\ 








i i i 


R 

c 

i 





-15 -10 -5 5 10 

ln(mr) 
FIG. 1. The electric field E in units of the critical field E c = m 2 /e is shown as a function of the 

distance from the string r for A = 4 x 10 3 m and 5 = 10~ 5 m -1 (solid line). Dotted lines indicate the 

analytic approximations ([!]) and ( Jl8| ) in the appropriate regimes. The core radius 5, the screening 

radius R s and the condensate radius R c are also indicated. 
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FIG. 2. The condensate radius R c vs the core radius 5 for A = 270m 1 . 
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FIG. 3. The condensate radius R c vs the linear charge density of the string, A for 5 = 0.001m -1 . 
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FIG. 4. The screening radius R s vs the linear charge density of the string A for 5 = 0.001m 1 . 
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